Resonance width distribution for open quantum systems 
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Recent measurements of resonance widths for low-energy neutron scattering off heavy nuclei show 
large deviations from the Porter- Thomas distribution. We propose a "standard" width distribution 
based on the random matrix theory for a chaotic quantum system with a single open decay channel. 
Two methods of derivation lead to a single analytical expression that recovers, in the limit of very 
weak continuum coupling, the Porter-Thomas distribution. The parameter defining the result is the 
ratio of typical widths T to the energy level spacing D. Compared to the Porter-Thomas distribution, 
the new distribution suppresses small widths and increases the probabilities of larger widths. We 
show also that it is necessary to take into account the gamma channels. 
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I. INTRODUCTION 

Random matrix theory as a statistical approach for 
exploring properties of complex quantum systems was 
pioneered by Wigner and Dyson half a century ago [l|. 
This theory was successfully applied to excited states of 
complex nuclei and other mesoscopic systems [2- 5] , eval- 
uating statistical fluctuations and correlations of energy 
levels and corresponding wave functions supposedly of 
"chaotic" nature. 

The standard random matrix approach based on the 
Gaussian Orthogonal Ensemble (GOE) for systems with 
time-reversal invariance, and on the Gaussian Unitary 
Ensemble (GUE) if this invariance is violated, was formu- 
lated originally for closed systems with no coupling to the 
outside world. Although the practical studies of complex 
nuclei, atoms, disordered solids, or microwave cavities 
always require the use of reactions produced by external 
sources, the typical assumption was that such a probe 
at the resonance is sensitive to the specific components 
of the exceedingly complicated intrinsic wave function, 
one for each open reaction channel, and the resonance 
widths are measuring the weights of these components 
Q. With the Gaussian distribution of independent am- 
plitudes in a chaotic intrinsic wave function, the widths 
under this assumption are proportional to the squares of 
the amplitudes and as such can be described, for v in- 
dependent open channels, by the chi-square distribution 
with v degrees of freedom. For low-energy elastic scatter- 
ing of neutrons off heavy nuclei, where the interactions 
can be considered time-reversal invariant, one expects 
v = 1 that is usually called the Porter-Thomas distribu- 
tion (PTD) 0. 

Recent measurements [H, [t| claimed that the neutron 
width distributions in low-energy neutron resonances on 
certain heavy nuclei are different from the PTD. As a 
rule, the fraction of greater widths is increased, while 
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FIG. 1: The proposed resonance width distribution accord- 
ing to eq. @ with a single neutron channel in the practically 
important case rj ^> V. The width T and mean level spacing 
D are measured in units of the mean value (r) . 



the fraction of narrow resonances is reduced which, be- 
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FIG. 2: The proposed resonance width distribution accord- 
ing to eq. (|27[) with a single neutron channel and N gamma- 
channels in the practically important case rj S> T. The neu- 
tron width r, radiation width 7, and mean level spacing D 
are measured in units of the mean value (r). 
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ing approximately presented with the aid of the same 
standard class of functions, would require v ^ 1. The 
literature discussing the scattering and decay processes 
in chaotic systems, see for example (Iol - fl3j and refer- 
ences therein, does not provide a detailed description of 
the width distribution for the region of relatively small 
widths as observed in low-energy neutron resonances. 

There are various reasons for possible deviations from 
the simple statistical predictions [T34l6l |. First of all, 
the intrinsic dynamics, even in heavy nuclei, can be dif- 
ferent from that in the GOE limit of many-body quan- 
tum chaos. If so, the detailed analysis of specific nu- 
clei is required. As an example we can mention 232 Th, 
where for a long time a sign problem exists [I?} con- 
cerning the resonances with strong enhancement of par- 
ity non-conservation in scattering of longitudinally polar- 
ized neutrons. The observed predominance of a certain 
sign of parity violating asymmetry contradicts to the sta- 
tistical mechanism of the effect and may be related to 
the non-random coupling between quadrupole and oc- 
tupole degrees of freedom [l8| . The width distribution in 
the same nucleus reveals noticeable deviations from the 
PTD. The presence of a shell- model single-particle res- 
onance serving as a doorway to the compound nucleus 
can also make its footprint distorting the statistical pat- 
tern. Another (maybe related to the doorway resonance) 
effect can come from the changed energy dependence of 
the widths that is usually assumed to be proportional to 
£ or neu -t;rons with orbital momentum I. Finally, 
the situation is not strictly one-channel, since, along 
with elastic neutron scattering, many gamma-channels 
are open as well. However, apart from structural effects, 
even in one-channel approximation, there exists a generic 
cause for the deviations from the PTD, since the appli- 
cability of the GOE is anyway violated by the open char- 
acter of the system . The appropriate modification of 
the GOE and PTD predictions, which should be applied 
before making specific conclusions, is our goal below. 

The resonances are not the eigenstates of a Hermitian 
Hamiltonian, they are poles of the scattering matrix in 
the complex plane. Their complex energies £ = E — iT/2 
can be rigorously described as eigenvalues of the effec- 
tive non-Hcrmitian Hamiltonian [201 ] - As shown long ago, 
even for a single open channel, the statistical properties 
of the complex energies cannot be described by the GOE. 
The new dynamics is related to the interaction of intrinsic 
states through continuum. In the limit of strong coupling 
this leads to the overlapping resonances, Ericson fluctua- 
tions of cross sections, and sharp redistribution of widths 
similar to the phenomenon of super-radiance, see the re- 
view 21( and references therein. The control parameter 



of such restructuring is the ratio 



k = 



ttT 
2D 



(1) 



of typical widths, T, to the mean spacing between the 
resonances, D. In the region of low-energy neutron reso- 
nances, n is still small but in order to correctly separate 



the general statistical effects from peculiar properties of 
individual nuclei we need to have at our disposal a generic 
width distribution that differs from the PTD as a func- 
tion of the degree of openness. 



II. RESONANCE WIDTH DISTRIBUTION 

The goal of the paper is to provide a practical tool 
that would allow one to compare an experimental out- 
put for an unstable quantum system with predictions 
of random matrix theory. We propose a new distribu- 
tion function that is based, similar to the GOE, on the 
chaotic character of time-reversal invariant internal dy- 
namics and corresponding decay amplitudes, but prop- 
erly accounts for the continuum coupling through the ef- 
fective non-Hermitian Hamiltonian. The numerical simu- 
lations for this Hamiltonian were described earlier [l5l[22l | 
but here we derive the analytical expression. We limit 
ourselves here by the situation typical for nuclear appli- 
cations, namely k < 1. The super-radiant regime, k > 1, 
can be of special interest, including such systems as mi- 
crowave cavities, and in the considered framework the 
formal symmetry exists, k — > 1/k. At a large number 
of resonances and fixed number of open channels, after 
the super-radiant transition the broad state becomes a 
part of the background while the remaining "trapped" 
states return into the non-overlap regime. However, in 
heavy nuclei this transition hardly can be observed be- 
cause earlier many new channels can be opened; in the 
modification of the PTD we see only precursors of this 
transition. 

Our arguments will follow two different routes which 
lead to the equivalent results. The final formula for the 
statistical width distribution can be presented as 



P(T) = C 



exp[-^r(^-r)] /smh 
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Here we consider iV 3> 1 intrinsic states coupled to a 
single decay channel, for example, s-wave elastic neutron 
scattering. The parameter D is a mean energy spacing 
between the resonances, k is a new dimensionless combi- 
nation, cq. ([1]), and C is a normalization constant. The 
quantity rj is the total sum (the trace of the imaginary 
part of the effective non-Hermitian Hamiltonian that re- 
mains invariant in the transition to the biorthogonal set 
of its cigenfunctions) of all N widths; it appears as a 
parameter that fixes the starting ensemble distribution, 
see eqs. ([8]) and ©. The possible values of widths are 
restricted from both sides, < T < rj. The above men- 
tioned symmetry k — > 1/k is reflected in the symmetry 
r — » rj — r of a factor in eq. @ but, as was already 
stated, our region of interest is at T <C r\. Another pa- 
rameter, (j, determines the standard deviation of variable 
r evaluated consistently with the distribution of eq. ((2|) . 
In the practical region far away from the super-radiance 
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we obtain 



p(r) = c Xl (r; 



I sinh i 



1/2 



(3) 



The PTD is recovered in the limiting case /t< 1 that cor- 
responds to the approximation of an open quantum sys- 
tem by a closed one. The new element is the factor explic- 
itly determined by the coupling strength n. With grow- 
ing continuum coupling the probability of larger widths 
increases. The distribution ((2J for different ratios (T) /D 
is shown in Fig. 1. 

The origin of the square root in the new factor is 
the linear energy level repulsion typical for the GOE 
spectral statistics. Indeed, in the complex plane, £ = 
E — iT/2, the distance between two poles £ m and £ n is 
\J{E m — E n ) 2 + (r m — T„) 2 /4; after integration over all 
variables of other states we obtain a characteristic square 
root in the level repulsion, see below eq. (fT2]) and the dis- 
cussion after eq. ([2"T]) . 



III. EFFECTIVE NON-HERMITIAN 
HAMILTONIAN AND SCATTERING MATRIX 

In order to come to the result ([2]), we start with the 
general description of complex energies £ = E — iT/2 in 
a system of N unstable states satisfying the GOE statis- 
tics inside the system and interacting with the single open 
channel through Gaussian random amplitudes. The gen- 
eral reaction theory [23| is constructed in terms of the 
elements of the scattering matrix in the space of open 
channels a, &, 



S ba (E) = S ba - iT ba {E). 



(4) 



Within the formalism of the 
Hamiltonian H = H - (i/2)W, 



effective non-Hermitian 
the T-matrix is defined 



T ba {E) 



N 

E 
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(5) 



in terms of the amplitudes A® connecting an internal 
basis state n with an open channel a. Here we neglect 
the potential part of scattering that is not related to the 
internal dynamics of the compound nucleus. The anti- 
Hcrmitian part of the effective Hamiltonian is exactly 
represented by the sum over k open channels, 



W„ 



a=l 



A a A a * 



(6) 



where the amplitudes can be considered real in the case 
of time-reversal invariance. It is important that the fac- 
torized structure of the effective Hamiltonian guarantees 
the unitarity of the scattering matrix. The amplitudes 



A°^ are uncorrelated Gaussian quantities with zero mean 

and variance defined as A a n A h n , = S ab S nn ^/N. The trace 
of the anti-Hermitian part of the effective Hamiltonian, 
i] = TrW, i.e. the total sum of all N widths used in 
eq. ([2]), is a quantity invariant under orthogonal trans- 
formation of the intrinsic basis. The detailed discussion 
of the whole approach, numerous applications and rele- 
vant references can be found in the recent review article 

The simplest version of the i?-matrix description uses 
instead of the amplitude T ba its approximate form, where 
the denominator contains poles on the real energy axis 
corresponding to the eigenvalues of the Hermitian part H 
of the effective Hamiltonian. Then the continuum cou- 
pling occurs only at the entrance and exit points of the 
process while the influence of this coupling on the in- 
trinsic dynamics of the compound nucleus is neglected 
(in general, H should also be renormalized by the off- 
shell contributions from the presence of the decay chan- 
nels). Contrary to that, the full amplitude T ba , eq. (|5|), 
accounts for this coupling during the entire process in- 
cluding the virtual excursions to the continuum and back 
from intrinsic states. The poles are the eigenvalues of the 
full effective Hamiltonian in the lower half of the complex 
energy plane. The experimental treatment corresponds 
to this full picture. According to the original paper 
the i?-matrix code SAMMY [24J had been used in the ex- 
perimental analysis where the relevant expression is given 
in the form 



7Ac7Ac' 



E - iTx/2 



Kjj> 



(7) 



and the treatment included a careful segregation of s- 
and p-resonances, J = J' = 1/2 for an even target nu- 
cleus. In the notations of [24| A represents a particular 
resonance, E\ is the energy of the resonance. Here we 
can identify the intermediate states A and their complex 
energies E\ — iT\/2 with the eigenstates and complex 
eigenvalues of "H, while the numerator includes the ampli- 
tudes transformed to this new basis (under time-reversal 
invariance the scattering matrix is symmetric). In terms 
of the reduced width 7^ c and the penetration factor P c , 
the partial width is T\ c = 2P c jf c . Assuming a single 
channel and universal energy dependence of penetration 
factors, the statistics of the total widths is the same as 
that of 7^ c . 



IV. FROM ENSEMBLE DISTRIBUTION TO 
SINGLE WIDTH DISTRIBUTION 

For a single-channel case, the joint distribution 
P(E; r) of all complex energy poles has been rigorously 
derived in [l9| under assumptions of the GOE intrinsic 
dynamics in the closed system and Gaussian distributed 
random decay amplitudes. The result is given by 
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P(E;f)=C N ] 



(-Br, 



E„ 



(E m — E n ) 



(r m +r„)= 

4 



n 



,-NF(E;T) 



(8) 



where the "free energy" F contains interactions of TV ^> 1 
complex poles in the interval 2a = ND of energies, 

F & r) = ^ E E « + i E r - r « + i E r -( 9 ) 

n m<in n 

For given AT, this distribution contains two parameters, 
the semicircle radius a for the intrinsic dynamics and n 
characterizing the total trace of the imaginary part of the 
effective Hamiltonian. 

Considering this free energy in the "mean-field" ap- 
proximation, we see that the original mean value (F)o = 
rj/N is substituted by (r) that is determined by the com- 
petition of two terms, l/(r) = 1/(T) + (r)/4D 2 . The 
first product in front of exp(— NF) in eq. (HJ substitutes 
the GOE level repulsion by the repulsion in the complex 
plane and interaction of the poles with their negative- 



r "images" [l9|. The structure of this result guarantees 
that all widths V are positive. The difficulty with the dis- 
tribution of eq. ((8]) is that it is not an analytic function 
of complex energies. 

Our first step is to specify a single 7V-th pole 
(En,Fn) = {E,T) and, using the fact that the distri- 
bution ensures T n > 0, return to the absolute values 
of the amplitudes, y/T^ — for other roots. In this 
form we can apply the steepest descent method owing 
to a large parameter AT 3> 1 and a saddle point inside 
the integration interval that was absent in the initial ex- 
pression. Integration Il^Li 1 d£n over ai l but one variable 
r = Tn = Cat leads to the following result for P(E;T) 
as a function of multiple energy variables E and a single 
width variable T: 



P{E; T) = C N ][ \E m - E n \ \[ J (E n - Ef + — exp 



r 2 



N 



exp 



■2V 



2tt 



N-l 



(10) 
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Introducing new variables, 2a /N = D and A = r)N, we factors in cq. (fTU)) in the limit of A^ —> 00: 
shall examine the behavior of one of the A^-depcndent 
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(11) 



The scaling properties of the parameters 77, a 2 , and A 
are as follows: n oc N, a 2 and A are both oc A^ 2 . As 
a result, the product of two exponents produces a well 
defined limit that brings in the desired dependence on 
the coupling strength k. 

The real energy distribution does not change much in 
an open system with a single decay channel being still, at 
finite but large N, close to a semicircle. We are working 



in the central region of the spectrum where the level den- 
sity is approximately constant and the energy spectrum 
is close to equidistant (the maximum of the level spacing 
distribution is always at s = SE/D w 1 although the dis- 
tribution in an open system changes at small spacings, 
s < 1, as we will comment later). With E n = E + nD, we 
are able to perform an exact calculation of the product: 



C N ]jJ(E n -EY 



£2 

T 



c 



N 



f N 

n 

\n=l 



r 2 /4 

{nD) 2 



1/2 
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V 2 D 
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(12) 
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where we have used the famous Eulcr formula, 



sinh x 



n 

fc=i 



(13) 



The width-independent factors will enter the normaliza- 
tion constant. Of course, the whole reasoning is valid 



in the limit N ^> 1. Finally, the width distribution for 
r <C r\ is represented by 



J 



/ sinh [f£ \ V2 ex P 

\ H ) Vr 



P(T) = C 



_/v-p 



exp 



2a 2 



(14) 



DOORWAY APPROACH 



of the secular equation, 



As an alternative derivation, we will apply the doorway 
approach 0, HH, [26[ • Here we use the eigenbasis of the 
imaginary part W of the effective non-Hermitian Hamil- 
tonian. Due to the factorized nature of W dictated by 
unitarity (l9j . the number of its non-zero eigenvalues is 
equal to the number of open channels. In our case we 
have only one eigenvalue, the doorway eg — irj/2, that 
has a non-zero width equal to the imaginary part 77 of 
the trace of the Hamiltonian. Remaining basis states are 
stable being driven by the Hcrmitian intrinsic Hamilto- 
nian; its diagonalization produces their real energies e n . 
These states acquire the widths through the interaction 
with the doorway state; the corresponding matrix ele- 
ments will be denoted h n . In this basis, the Hamiltonian 
is represented as 



(15) 
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The complex eigenvalues 8 = E — ir /2 are the roots 



i ^ \h n \ 2 



(16) 



that is equivalent to the set of coupled equations for real 
and imaginary parts, 



N 



E = e + J2\ h n\ 



E — e n 



(£- £n ) 2 + r 2 /4' 



(17) 



r = 



n 



= f(T,E). (18) 



(E-e n ) 2 +r 2 /4 



For the Gaussian distribution of the coupling ma- 
trix elements with ( \h\ 2 ) = 2<r 2 /N (this scaling was 
derived in HH), we obtain 



P(T) 



+ OO 



(5(r-/(r,£))exp 



cr 2 ^ r ' 
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N 



I dh r , 



(19) 



r 



The integration in ([19]) via the steepest descent method 
leads to eq. $Z§. In order to get this result we use a pos- 
sibility to find a highest root E = ejy which we set as an 
origin relative to which the energies e n can be counted as 



E — En, £n — £n —nD. An important intermediate step 
is the evaluation of the infinite product of the Lorentzian 
peaks that can be simplified as 
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--N) 



-1/2 
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-1/2 
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(20) 



In a similar way one can analyze the resonance spac- 
ing distribution P(s) along the real energy axis; spacings 
s = 6E/D are measured in units of their mean value 
D. As predicted in (l9j and observed numerically in [22| . 
the short-range repulsion disappears and the Wigner sur- 
mise with the standard linear preexponential factor s is 
substituted by the square root, 



p(a) 



4Jcxp[- 



const • s 2 ] . (21) 



At spacing s«l, the probability behaves as a + bs 2 with 
the quadratic dependence on s that, similar to the GUE, 
mimics the violation of time-reversal invariance due to 
the open decay channel. The absence of short-range re- 
pulsion, a (the interaction through continuum, oppo- 
site to a normal Hermitian perturbation, repells widths 
and attracts real energies J27| ) , reflects the energy uncer- 
tainty of unstable states. 

We demonstrated that two complementary approaches 
which reflect different physical aspects of the situation 
lead essentially to the equivalent (after identification of 
corresponding parameters) results which we prefer to 
write in the form (|2J). We expect that for other canoni- 
cal ensembles the width distribution far from the super- 
radiance can be expressed by a similar formula with the 
function (sinh k/k)' 3 / 2 , where the standard index of en- 
semble is P = 1 for the GOE, (3 = 2 for the GUE, and 
j3 = 4 for the Gaussian Symplectic Ensemble. In the 
same way we expect the square root in eq. (|2"Tj) to be 
substituted by the same power (3/2. 

The doorway approach naturally indicates the limits 
of the variable, < T < i). It has also an advan- 
tage of the possibility to generalize the answer taking 
into account explicitly the rigidity of the internal energy 
spectrum with fluctuations of level spacings around their 
mean value D [in our approximation, only this average 
value enters eq. (0)]. Another direction of generaliza- 
tion includes the possible influence of a single-particle 
resonance depending on a position of its centroid with 
respect to the considered interval of the resonance spec- 
trum. In particular, that centroid may be located under 
threshold of our decay channel. In this case even the 
standard energy dependence of the widths can change as 
was mentioned long ago [25j], see also [l4j]. The doorway 
state may or may not coincide with such a resonance so 
that the effective Hamiltonian (|T5|) may contain two spe- 
cial states coupled with the "chaotic" background, one by 



intrinsic interactions and another one through the con- 
tinuum. 

VI. ADDING GAMMA-CHANNELS 

The goal of this section is to estimate in the same spirit 
the influence of 7-channels on the resonance width distri- 
bution. Only a single open elastic neutron channel was 
taken into account in the analysis of data [H, |{|. The 
presence of even weak additional open channels changes 
the unitarity conditions. Examples of mutual influence 
of neutron and gamma channels are well known in the 
literature from long ago, see for example [28| . 

Generalizing the doorway description we allow now 
each intrinsic state to decay by gamma-emission which 
is always possible independently of the position of the 
neutron threshold. In the simplest approximation, the 
effective non-Hcrmitian Hamiltonian is now represented 
by 
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(22) 



where we assumed that the intrinsic part of the matrix 
is pre-diagonalized and introduced 7„ as the widths for 
7-channels. 

Analogously to eq. (frB]) . the complex energy eigenval- 
ues £ = E — are the roots of the secular equation, 
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\h„\ 



~1 & { £ n 2 ^ n ) 

equivalent to a set of coupled equations, 

E — e n 



(23) 
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E = e + ^2\h r , 



n=l 



(E - e n y + (r - 7n )2/4 



. (24) 



i + 2^„=ll'M (_E-e„)^ + (r- 7 „)V4 

The resonance width distribution for an open quantum 
system with 7-channels included is given by 
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Estimating the gamma-widths by their average value, single open channel we come to the final expression for 
7, and acting in the same manner as in the case of a the resonance width distribution, 

I 
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(27) 



that is shifted by T — > T — 7 compared to the previ- 
ous result. The mentioned earlier symmetry between the 
ends of the distribution, r = and r = 77, would be 
substituted here by T — >• (rj + 7) — T. Thus, the effective 
influence of 7-channels on the resonance width distribu- 



J 



tion is reduced here to a shift of the whole distribution 
by a mean radiation width 7 as seen in Fig. 2. In the 
practical region far away from the super-radiance, T <C 77, 
we obtain 



p(r,7) = x ?[r- 7 ] 



sinh 


'f(r-7)" 


2D 


TV 


T-7) 


2D 



1/2 



(28) 



In order to extract the neutron width from the total reso- 
nance width, the treatment of the data has to be modified 
making in a sense an inverse shift. Of course, a more pre- 
cise consideration should use a statistical distribution of 
the gamma widths. 



"super-radiant" state becomes essentially the part of the 
background, while the remaining "trapped" states return 
to the weak coupling regime. This physics was repeatedly 
discussed previously, especially in relation to quantum 
signal transmission through mesoscopic devices (2lL [2{| , 
but it is outside of our interest here. 



VII. CONCLUSION 

In this article we propose a new resonance width dis- 
tribution for an open quantum system based on chaotic 
intrinsic dynamics and coupling of states with the same 
quantum numbers to the common decay channel. Two 
approximate methods lead to an equivalent analytical ex- 
pression for the width distribution that does not belong 
to the class of chi-squarc distributions with the only pa- 
rameter v traditionally used in the analysis of data. In 
the limit of vanishing openness and return to a closed 
system we recover the standard PTD. The new result 
depends on the ratio (pj of the width to the mean level 
spacing, k ~ T/D, that regulates the strength of the 
continuum coupling. The deviations from the PTD grow 
with k up to the critical strength k ~ 1, when the broad 



In the practical region of low-energy neutron reso- 
nances, the effects predicted here are relatively small. 
Although at small k the derived neutron width distribu- 
tion supports an experimental trend, the final judgment 
can be made only after the presence of gamma-channels 
was accounted for. We have to attract the attention of 
experimentalists to the fact that the data should be an- 
alyzed with the aid of the distribution that does not be- 
long to the routinely used chi-square class; gamma chan- 
nels should be included into consideration. We can also 
mention that the result agrees with numerical simula- 
tions [lH for the full many-resonance distribution func- 
tion ([5]). Using the suggested distribution as a new refer- 
ence point, one can ascribe the remaining deviations to 
the specific features of individual systems (level densities, 
single-particle structure in a given energy region, shape 
transformations, energy dependence of the widths etc.). 
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Unfortunately, we still do not have experimental tests for 
the full distribution ©. Although in nuclear physics it 
is hard to make such a detailed analysis for higher ener- 
gies and greater degree of resonance overlap, the systems 
with tunable chaos, such as microwave cavities, acous- 
tic blocks, or even elastomechanical devices [3(|, seem to 
provide appropriate tools for such studies. 
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